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LM Justification du choix du suiet 


Les réseaux d'assainissement au Maroc, en Europe en Angleterre et même aux Etats Unis 
connaissent aujourd'hui une situation de dégradation. Les manifestations de cette situation, 
variant selon les pays, apparaissent à travers le mauvais état des équipements, leurs faibles 
performance et efficacité. Mais le problème de la dégradation va au-delà de ces 
manifestations de nature physique. Il réside dans le fait que les capacités réelles des 
équipements ne s'avèrent plus compatibles avec les exigences de qualité et les objectifs du 
service. 


Les équipements existants sont souvent trop petits par rapport à cet objectif de traitement 
et rendent nécessaire la construction de nouvelles infrastructures. Si de nouveaux 
investissements sont nécessaires, ils doivent être faits dans des conditions qui semblent 
économiquement acceptables pour les autorités locales. Pour limiter le nombre de 


nouvelles infrastructures, on cherche à gérer de façon optimale les capacités de traitement 
de pollution par une meilleure utilisation des infrastructures existantes (collecteur, bassin 
de rétention, déversoir). 


Il est donc clair que les réseaux d'assainissement apparaissent comme des systèmes 
complexes de tracés et de matériaux dont l'optimisation ne peut être réalisée par une 
démarche énumérative. 


La lutte contre les inondations demeure un problème réel dans la majorité des grandes villes 
du monde. Si l'eau est source de vie, elle peut être également source de nuisances. En effet la 
forte urbanisation caractérisée par la minéralisation des espaces, la gestion de 
l'assainissement urbain a pris de plus en plus d'importance majeure depuis plusieurs années. 
Cette forte urbanisation a pour conséquence l’augmentation du ruissellement urbain vers les 
réseaux d’égout. La garantie d’un environnement sain et sécurisé aux riverains, devient un 
défi majeur avec la volonté d’évacuer le volume ruisselé le plus rapidement possible. A ces 
eaux pluviales s'ajoutent les eaux usées produites par les populations elles-mêmes qu'il 
faudra aussi rapidement évacuer et traiter quotidiennement dans les normes 


I. Problématique 
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Dans l'étude des réseaux d'assainissement /a première étape peut, par exemple, consister à 
vérifier qu'une manipulation des organes de contrôle tels que les vannes et pompes est 
capable de minimiser les déversements vers le milieu naturel. Cette gestion, que l'on 


appellera « gestion de référence », permet de déterminer les stratégies de commande sur 


toute la durée de l'événement pluvieux connu à l'avance. Ce calcul se fait donc à la fin de 
l'événement pluvieux et permet de dire ce qui aurait pu être fait avec les organes de régulation 


en terme de minimisation des volumes déversés. par les graphes 





et la programmation linéaire mixte permettent de déterminer une solution optimale. 
Cet article s'intéresse à la vérification de l'optimalité et à l'applicabilité de la programmation 





linéaire par les graphes comparée à la programmation linéaire mixte dans le cas de la « gestion 
de référence » sur le réseau d'assainissement (Saverne (France) comme exemple). 


I. Objectif du sujet 


L'objectif de cet article est, dans un premier temps, de vérifier l'optimalité de la 


programmation linéaire par les graphes dans laquelle on impose au graphe un ordre de 
saturation des arcs pour prendre en compte l'hydraulique du réseau que l'on détaillera par la 


suite. Cette vérification se fait par rapport à la programmation linéaire mixte utilisant un 
algorithme basé sur la « branch-and bound method » capable de prendre en compte les 
variables binaires traduisant les contraintes de saturation précédentes. Pour les deux 
approches, on utilisera la même formulation mathématique. La comparaison se fera dans le 
cas de la « gestion de référence » sur le réseau d'assainissement de Saverne (France) et 
porte sur les volumes déversés vers le milieu naturel. Dans un deuxième temps, on compare 
les temps de calcul de chaque méthode. On pourra ainsi vérifier l'applicabilité des 
techniques de calcul en vue de la gestion en temps réel d'un réseau d'assainissement. Le but 
final est de pouvoir choisir un compromis entre la pertinence ou la précision des résultats 
et le moyen de les obtenir. Nous allons, dans un premier temps, décrire le principe de 
modélisation hydraulique des principaux ouvrages (collecteur, bassin d'orage et déversoir 
d'orage) par la programmation linéaire par les graphes et par la programmation linéaire 
mixte. Dans un deuxième temps, après une présentation succincte du réseau 
d'assainissement de Saverne, nous détaillerons les résultats des optimisations. 


v. Modélisation Mathématique des ouvrages 
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1. Système d'équations à résoudre 


Selon la programmation linéaire par les graphes ou mixte, on peut écrire le 
sous la forme : 


Min (2 5-1 ii) 





Sous les contraintes : M. @ = B et 0 < p < cm. 


n: le nombre d'arcs (Arc : transfert d'un volume entre deux lieux du réseau 


d'assainissement), Yi : coût unitaire de transit, @;: flot (volume) de transit, 


M : matrice d'incidence, B vecteur dépendant du site d'étude et cm capacité 
de transit maximale. Le calcul consiste en la minimisation d'une fonction linéaire 
sous un certain nombre de contraintes linéaires avec variables réelles et binaires. La 
fonction linéaire représente le volume déversé minimisé sur l'ensemble de la durée 
de la pluie. Un coût unitaire (Constant) par volume n'est affecté qu'aux arcs 
déversants vers le milieu naturel. Les variables réelles sont les volumes d'eau écoulés 
dans chaque collecteur pendant un pas de temps donné. Les contraintes et les 
variables binaires représentent les lois hydrauliques des réseaux d'assainissement 
que nous détaillerons par la suite. On en déduit, par exemple, les temps de 
démarrage/d'arrêt de pompes ou les temps d'ouverture/fermeture de vannes. 


Le pas de temps doit être assez court pour que la gestion soit réalisée en temps réel 
et soit suffisamment précise pour faire face aux pluies violentes de courte durée. Il 
doit être aussi assez long pour obtenir toutes les informations des différents capteurs 
de débit, par exemple (mesure de l'état du réseau), pour calculer les nouvelles 
consignes et pour les envoyer à chaque ouvrage de contrôle (pompes, vannes, ...). 
Pour déterminer ce pas de temps, il faut aussi connaître le temps de transfert des 


débits dans le réseau. Pour ces raisons, le pas de temps choisi sur le réseau de 


Saverne est de 4 minutes. 


2. Propriété d'implication 
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Nous aurons besoin de traduire en programmation linéaire mixte les 
implications 








que nous verrons par la suite. I| faut écrire sous forme de 
contraintes linéaires des conditions de saturation du type : « si f(s) > 0 alors 


g(s) = gmax ». L'idée est ici de trouver des fonctions f et g qui respectent les 
conditions d'écoulement. 


Condition : f(s) > 0 => g(s) = gmax avec pour contrainte 


0 < f(s) < fmax et 0 < g(s) < gmax Pour transformer cette implication en 
contraintes qui pourront alors s'ajouter à celles du programme linéaire, nous 
introduisons une variable binaire supplémentaire x. On écrit alors : 


f(s) < fmax (1 - x) et g(s) > gmax (1 - x). Nous obtenons bien l'implication 
recherchée puisque si f(s) > 0, x = 0 (la valeur 1 étant impossible puisqu'elle 
donne f(s) < 0) et on a g(s) > gmax et compte tenu de la contrainte g(s) < gmax 
Alors g(s) = gmax. Il faut aussi remarquer que, quand f(s) < 0, x peut prendre la 
valeur 0 ou 1. Par contre si g(s) < gmax alors x = 1 et dans ce cas 


f(s) = 0. Nous avons utilisé cette technique pour modéliser l'écoulement à 
travers les ouvrages. 


3. Modélisation des ouvrages 


Dans cette partie, nous allons décrire la modélisation du comportement hydraulique des 
collecteurs, déversoirs d'orage et bassins d'orage. Nous donnerons également les limites de 
chaque modèle. 


3.1 Modélisation de l'écoulement dans un collecteur 


L'écoulement dans un collecteur est représenté mathématiquement par les équations 
différentielles non-linéaires de Barré de Saint Venant en unidimensionnel. Compte tenu de la 
difficulté de résolution de ces équations, nous avons préféré utiliser le modèle de 
Muskingum (CUNGE, 1969) qui en est une simplification. 
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Certains phénomènes hydrauliques tels que les mises en charge ou les influences aval 
provoquées par exemple par la présence d'un seuil à l'aval ne seront pas pris en compte. 
Toutefois lorsque ces phénomènes hydrauliques ne sont pas importants, ils jouent un rôle 
qui peut être négligeable. On se place donc dans ces conditions hydrauliques. Le modèle de 
Muskingum utilisant les équations de stockage et de conservation propose une 
relation du type : 


Qs(t + At) = 6..Qelt) +62. Qe(t + At) +6; . Qs(t) où, Qs(t), Qe(t) représentent 
respectivement le débit de sortie et d'entrée à l'instant t, et 6,62, B3 des 
paramètres définis comme suit : 


2.K.a.+At —2.K.a.+At __2.K(1-a)-At 


P1 = 2.K.(1 — à) + At B2 = 2K(-@+At P3 = 2.K.(1 — à) + At MECS 


B1+B2+63 = 1, K et a sont des paramètres de calage du modèle de 
Muskingum modifiés (VASQUEZ, 1999). Afin de représenter le fonctionnement 
hydraulique d'un collecteur par la programmation linéaire par les graphes, on 
fait une décomposition spatiale (entrée/sortie du collecteur) et temporelle 
(décomposition en pas de temps) comme indiquées sur la figure 1. 


Figure 1 à gauche : Représentation par les graphes d'un transfert et d'une 
diffusion de débit. 


Figure 2 à droite : Décomposition d'un hydrogramme. 


A 
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Entrée du collecteur Sortie du collecteur Volume d'entrée Volume de sortie par ps 
ti G 0 t par pas de temps de temps 
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0 b 
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t 0 ot JA, M 
ee bou bthtbtktkt 
Décomposition spatiale 


L'arc 13 (Arc13) montre qu'un volume arrivant à l'instant t, se propage avec un temps de 
transfert de t1-t3. On voit également sur la figure 1 qu'un volume arrivant à l'instant t1 va 
se diffuser à travers les arcs Arc13, Arc14 et Arc15. La représentation sous forme 
d'hydrogramme de la figure 1 peut être représentée par la figure 2. Le volume V1 arrivant à 
l'instant t1 va être transféré par les volumes dans les Arc13, Arc14 et Arc15. On fait de 
même pour les volumes V2 et V3. 


Pour déterminer le volume des arcs, on utilise la technique de l'hydrogramme unitaire. Par 
cette technique de décomposition spatiale et temporelle, on prend en compte la diffusion 
et le temps de transfert des débits en réseaux d'assainissement. On constate ainsi que dans 
un collecteur plus le débit est faible, plus le temps de transfert est long. Afin de prendre en 
compte ce principe, les arcs les plus éloignés sont d'abord saturés (c'est-à-dire qu'ils 
atteignent la capacité maximale en terme de volume) dans le sens 3-2-1 de la figure ci- 
dessous qui représente le Principe de saturation des arcs. 
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t 
to O 
t3 O a 
t4 O g> 







O Pour fixer la capacité maximale de 
chaque arc, on utilise le principe décrit par la figure 2. On simule par le modèle 
de Muskingum N hydrogrammes rectangulaires. On détermine la capacité 
maximale (volume d'eau maximal transité) pour chaque arc Arc;i;:2, ATCii+3 
et ArC;;44 partant de l'instant t;, de telle sorte qu'elle soit constante quel que 
soit l'instant tį. 


L'algorithme COLLECTEUR : Formulation pour la théorie des graphes : Les arcs 


doivent être saturés dans l'ordre décroissant des temps de transfert, d'abord 
i = n puis n-1 jusqu'à 1. Donc, saturer les n tronçons dans l'ordre (d'abordi=n, 
i=n-1i= 1) revient à dire d'un point de vue mathématique que : 


Si Ji€ {1, ..., n} tel qu'un volume transitant dans l’ arc; > 0 alors le volume 
transitant dans l’ arc;,,= capacité maximale de l’ arc;,..C'est-à-dire : 


1i€ {1, ..., n} tel que d;> 0 alorsd;,, = cm;,.. Avec di1 : le Volume transitant 
dans l'arc i + 1, et cm;., : la capacité de transit maximale dans l'arc i + 1. 
Formulation pour la programmation linéaire mixte : 


V i, 0 < d; < cm: 
di, 2 cm 1 +1 (1 — Xi, 1) x : Une variable binaire. 


3.2 Modélisation d'un bassin d'orage 
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La modélisation du volume d'eau dans un bassin d'orage se base sur la 
conservation de la matière : Vpassin () - Vpassin(t — 1) = Vsortant (t) 
+V'aéversé (t) - Ventrant (t) comme la montre la figure ci-dessous. 
Représentation élémentaire du bassin d'orage. Elementary représentation of 
the Storm water basin. 








V bassin (t 1 ) y Vies (t) 





néant (t) | Vaaran (t) 


V bassin (t) 


Ce schéma est reproduit à chaque instant, ce qui permet de fournir la 
décomposition temporelle du volume d'eau dans le bassin. Il n'existe qu'une 
seule implication à vérifier sur un bassin d'orage : le bassin ne peut déverser 
qu'une fois plein. 


L'algorithme BASSIN : Formulation pour la théorie des graphes : 


Si Ji€ {bassin} tel qu'un volume transitant dans l'arc, de déversement > 0 alors 
le volume transitant dans l’ arc;,.= capacité maximale du bassin. 


C'est-à-dire : Si Ji€ {bassin} tel quedbs;> 0 alors db;= cmb. Avec cmb : la 
capacité maximale des arcs correspond au volume total du bassin, db : le 
volume transitant dans l'arc i du bassin, dbs : le volume qui est rejeté vers le 
milieu naturel. 


Formulation pour la programmation linéaire mixte : 
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Vie {bassin? 0 < db; < cmb ; dbs; > 0 
dbs, < CP (1 - x) 
a > — 3 
db, 2 cmb (1 X;) CP : représente la 
capacité maximale de déversement pour un pas de temps. 


Déversement 


Coût non nul 


Milieu naturel 


Entrée V. Sortie 


Représentation du 
bassin d'orage par les graphes. Graphical representation of the storm water 
basin. 


3.3 Modélisation d'un déversoir d'orage 


La courbe de fonctionnement d'un déversoir peut être prise en compte par la 
programmation linéaire par les graphes en décomposant le processus de la 
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façon suivante : en appelant Q;éférence le débit à partir duquel le déversoir 
commence à déverser, la figure ci-dessous montre le principe de 
fonctionnement du déversoir. Les arcs sont saturés dans l'ordre suivant : 1, 2, 3 
et 4. 


Q 













Q 


déversé 


Q, ortie 


Q référence 


Qentrant 





C RAR 


4 Sortie déversée 
2 
Entrée EE, Sortie conservée 
3 
la courbe de fonctionnement du déversoir d'orage par les graphes. Graphical 
représentation of the storm overflow characteristics curve. 


Représentation de 


L'algorithme DÉVERSOIR : Formulation pour la théorie des graphes : 





Les arcs doivent être saturés dans l'ordre croissant de la numérotation des arcs 
du déversoir. Donc, saturer les n tronçons dans l'ordre (d'abord i = 1, i= 2, ...,i 
= n) revient à dire d'un point de vue mathématique que : 
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Si Ji€ {1, ...,n} tel qu'un volume transitant dans l’ arc; > 0 alors le volume 
transitant dans l’ arc;_, > 0 = capacité maximale de l'arc;_: du déversoir. 
C'est-à-dire : Si 4i€ {1, ..., n} alors 


d; > Otqad;_.=cmd;_.cmd;_,:lLa capacité maximale de l'arc i-1 
correspond au volume de déversement. 


Formulation pour la programmation linéaire mixte : 

vi O = dd, = cmdv. 

di < cmdv; (1-x;_.) 

d- = cmdv.., (1-x._.) 

V. Algorithmes de résolution utilisés 
L'algorithme « Branch and Bound » (séparation - évaluation) : 


L'intérêt de celui-ci est de s'affranchir d'un examen de toutes les solutions (qui 
représenteraient 2” possibilités où n est le nombre de variables binaires) tout en 
s'assurant d'obtenir la solution optimale. Afin de fournir une solution plus rapide sans 
pour autant aller jusqu'à la convergence exacte, le logiciel CPLEX encadre la fonction 
objective par des valeurs minimale et maximale pouvant être atteintes. La valeur 
minimale que l'on appellera « PLM min », consiste à fixer une partie des variables 
binaires comme étant des variables réelles comprises entre 0 et 1, alors que la valeur 
maximale que l'on appellera « PLM max », garantit les variables binaires comme étant 
égales à 0 ou 1 mais la solution trouvée est un minimum local. Cette démarche permet 
de stopper le calcul alors que le minimum global n'est pas atteint. En effet, dans 
certaines situations, On peut atteindre 24 heures de temps de calcul. L'encadrement du 
minimum global est donc assuré dans le cas où les temps de calcul deviennent très 
importants. 


L'algorithme primal-dual de Sakarovitch (Programmation linéaire par les 


graphes ) : calcul du flot maximal à coût minimal 
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Nous appellerons « PLG » la valeur de la fonction objective après convergence. Les modèles 
que nous avons développés précédemment imposent un ordre de saturation des arcs du 
graphe. Or, une des hypothèses à respecter pour garantir un minimum global en utilisant 
l'algorithme primal-dual est l'indépendance du choix des arcs. On ne garantit donc plus un 
minimum global. Les temps de calcul sont de l'ordre de 5 minutes. La partie suivante 
s'intéresse à la comparaison de la valeur « PLG » avec les valeurs de « PLM min » et « PLM 


max » pour les 34 pluies définies dans l'exemple de Saverne ci-dessous. 


vi. Apalication au réseau de Saverne 
Schémas du réseau : 


BASSIN D'ORAGE 


MONSWILLER 


HYDROGRAMME DO 





























3 MILIEU NATUREL 
MONSWILLER 
Groupe de 
pompage 
MILIEU NATUREL 
HYDROGRAMME y|régulation|; a 20 ma 
1 
SAVERNE Station de 
BASSIN D'ORAGE | traitement 
CHAMP DE FOIRE 
MILIEU NATUREL 
MILIEU NATUREL 
HYDROGRAMME pompage 
2 
SAVERNE SUD 


BASSIN D ORAGE 
ZORNHOFF 
Le Syndicat d'Assainissement de la Région de Saverne-Zorn-Mossel regroupe 10 communes 
situées à une distance de 35 km au Nord-Ouest de Strasbourg et compte environ 16 800 
habitants. Le réseau est principalement unitaire. || comporte 92 déversoirs d'orage, 12 
bassins d'orage et 3 stations de refoulement. Après simplification du réseau 
intercommunal, la figure en dessus de ce paragraphe présente le réseau d'assainissement 


Cou Liga) ds amod) du y Lol 


Ecole Mohammadia d'Ingénieurs 








= Lyl pull 344 el 
Université Mohammed V de Rabat 


OPTIMISATION DANS LES RESEAUX D'ASSAINISSEMENT (Gestion en temps réel) 
La théorie des graphes et la programmation linéaire mixte 
Préparé par : LAMINE Mohamed Encadré par : MME EL MKHALET 


étudié. Le bassin d'orage 4 (champ de foire) est régulé par l'intermédiaire d'une vanne 
motorisée alors que les bassins 11 (Monswiller) et 5 (Zornhoff) disposent chacun d'un 
groupe de pompage. 








volumes déversés vers le milieu naturel, nous imposerons donc un coût sur les volumes 
déversés. Le pas de temps choisi sur Saverne est de 4 minutes. Les événements pluvieux 
proviennent d'une base de données regroupant environ 2000 hyetogrammes. Les pluies 
ayant une lame d'eau totale inférieure à 2 mm provoquent un volume total inférieur à 4 % 
du volume du bassin ne sont pas prises en compte. Il reste cependant 535 événements 
pluvieux. Plusieurs paramètres ont été choisis pour extraire de cette base de données les 
événements pour l'optimisation de référence : 


e Événements provoquant un déversement dans le cas d'une gestion statique, 


e Durée de l'événement, + Apparition du débit de pointe en début, milieu et fin 
d'événement pluvieux, 


e Nombre de pics en débit par événement pluvieux, + Décalage temporel d'apparition des 
pointes de débit entre les 3 bassins versants. 


Compte tenu de ces critères, nous avons retenu 34 événements pluvieux. 


Résultats 
Différence Différence Différence Différence 
Pluie total des déversements des déversements des déversements des déversements 
Pluie + ui (mê) par rapport par rapport par rapport 
(m°) VoluMepig - au volume au volume au volume 
Volume, y max de la pluie des bassins total déversé 
1-50 13 031,1 40,4 0,30 % 0,96 % 6,28 % 
1-63 10 206,0 0,5 < 0,01 % 0,01 % < 0,01 % 
10-63 13 295,5 35,3 0,13 % 0,84 % 0,24 % 
10-85 16 589,6 5,7 0,02 % 0,14 % 0,03 % 
14-39 14 368,1 34,3 0,23 % 0,82 % 5,26 % 
15-25 9 099,3 13,1 0,12 % 0,31 % 0,79 % 


Écarts absolus et relatifs entre les déversements pour chaque bassin par événement pluvieux. 


-Les résultats du tableau 1 montrent que pour 6 pluies où la programmation linéaire mixte donne une 
meilleure solution, la différence des déversements peut atteindre 40 m3. Mais ceci ne représente 
que 0,30 % du volume total de la pluie, 0,96 % du volume total des 3 bassins et 6,28 % du volume 


déversé. 
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-Nous nous intéressons maintenant à la répartition des déversements selon les trois bassins. Pour le 
bassin 4, les écarts (Volume graphe -VOlUME jinéaire mixte) sont positifs ou nuls alors que pour les 2 


autres bassins, ils sont négatifs. Nous avons donc calculé le taux de déversement par bassin et pour 
chacune des pluies (tableau 2), c'est-à-dire pour la pluie p et le bassin n : 





% de déversement par % de déversement par 


pluie la programmation linéaire mixte la programmation linéaire par les graphes 
Bassin 4 Bassin 5 Bassin 11 Bassin 4 Bassin 5 Bassin 11 

1-50 100 % 100 % 

1-63 85 % 7 % 8 % 88 % 5 % 7 % 
10-63 91 % 0 % 9% 92 % 8 % 
10-85 84 % 9 % 7 % 84 % 9% 7 % 
14-39 100 % 100 % 

15-25 100 % 100 % 


Pourcentages de déversements par bassin et par événement pluvieux. 


On constate, comme le montre le tableau ci-dessus, que les déversements sont légèrement 
privilégiés à travers le bassin 4 pour la programmation linéaire par les graphes par rapport aux 
résultats obtenus par la programmation linéaire mixte. Ce constat s'est également révélé pour 
l'ensemble des autres pluies. Le fait que le bassin 4 soit toujours plus rempli et qu'il déverse toujours 
plus dans le cas de l'approche par les graphes peut venir de la numérotation. En effet, intéressons- 


nous au marquage des sommets et à la numérotation des arcs du graphe représentée à la figure 
suivante. 
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vers bassin d'orage 4 
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1000 
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1001 © 
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Fi 9 1500 | 
igure vers bassin d'orage 11 
Numérotation des arcs du graphe. 
O t= 500 





Arc numbering in the graph. 


Résultats Les 1500 premiers arcs ont pour origine le sommet 1 qui est l'entrée du système. 
Les 500 premiers arcs numérotés 1-500 sont les arcs d'entrée du bassin versant 1 (bassin 
d'orage 4), les arcs numérotés 501-1000 sont les arcs d'entrée du bassin versant 2 (bassin 
d'orage 5) et finalement les arcs numérotés 1001-1501 sont les arcs d'entrée du bassin 
versant 3 (bassin d'orage 11). L'approche par les graphes cherche d'abord à saturer le 
maximum d'arcs sans déversement. En effet, nous n'avons affecté qu'un coût aux arcs 
déversant vers le milieu naturel. Donc, l'algorithme va d'abord faire transiter l'eau vers la 
station de traitement (voir Schémas du réseau). Etant donnée la numérotation, elle traitera 
en premier l'eau du bassin versant 4. De même, les arcs de stockage du bassin 4 portent les 
numéros pairs de 23444 à 24440, ceux du bassin 5 les numéros pairs de 24444 à 25440 et 
ceux du bassin 11 les numéros pairs de 25444 à 26440. Une fois que, pour chaque bassin 
versant, le flot maximal à coût nul est calculé, c'est-à-dire qu'on ne peut pas envoyer plus 


d'eau à la station de traitement alors le flot supplémentaire est stocké. Les premiers arcs 


trouvés par l'algorithme sont fonction de la numérotation et correspondent au bassin 4. Il en 
est exactement de même pour les déversements. Ceci permet d'expliquer les constatations 
faites sur les déversements et les volumes stockés. 
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vi. Application au Maroc (cout, moyens...) ? C'est 
possible ? 


Comme vous le savez Les récentes inondations dramatiques de Taroudant, d'Er- 
Rachidia et de Khénifra ainsi que les différents commentaires et les appels à l'aide sur 
les réseaux sociaux ont déclenché la réaction du Gouvernement de Saäd Eddine El 
Othmani. 


A cet effet, le Ministère de l'Intérieur a élaboré un projet de lutte contre les 
inondations, dont Abdellah Nassif, Gouverneur dudit Projet a expliqué clairement que 
sur le plan technique, ce projet, conçu sur un système intégré, « recevra toutes les 
informations et les prévisions émises par les différents secteurs ». Ceux-ci se 
composent du département ministériel chargé de l'Eau et des agences de bassins 
hydrauliques, du secteur de l'urbanisme qui détermine les zones urbanisables, des 
collectivités territoriales ainsi que du Ministère de l'Intérieur qui agira en 
Coordinateur des opérations. Par ailleurs, le spécialiste indique que ledit système 

« permettra donc de mesurer les niveaux d'eau dans les villes, les villages et les 
douars », sans oublier de mentionner la réalisation d'une cartographie des 
installations et des habitants afin « de prévoir et gérer les risques 

d'inondations ». Ce projet, au regard des études, du système d'exploitation, des 
équipements électroniques, du centre au Ministère de l'Intérieur ainsi que des centres 
provinciaux, coûtera un montant de 32 millions de dirhams. 


De ce fait, il est clair que le Maroc aura besoin comme la plupart des autres pays 
d'une implantation des systèmes d'optimisation du type ‘optimal control”. En 
plus il a toutes les moyens nécessaires pour ce fait à savoir les grands 
ordinateurs, les équipements hydrauliques etc. 


VI. perspectives pour un élargissement de ce sujet 


Face au changement climatique, à une urbanisation croissante et une réglementation 


de plus en plus stricte, la gestion des réseaux d'assainissement s'avère toujours plus 
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complexe. Ainsi le développement de nouvelles solutions capables de contrôler les 
réseaux d’eau (pluviales+usées) en temps réel est nécessaire. Pour cela plusieurs 
logiciels très développés sont mis en disposition à savoir 

AQUADVANCED® Réseaux d'Eau qui permet de mesurer le flux, la pression, le 
débit et ainsi d'identifier et d'anticiper les anomalies intervenant sur le réseau d'eau 
potable. 

En plus de la collecte de ces informations, le logiciel assemble de très nombreuses 
données en provenance de sources variées comme le système de télégestion et 
d'acquisition de données SCADA (Supervisor Control And Data Aquisition), le SIG 


(Système d'Information Géographique) ou encore la télérelève. Il est devenu donc 


un véritable outil de pilotage des réseaux d’assainissement. 


AQUADVANCED® Assainissement permet de mobiliser les capacités de stockage au 
bon endroit et au bon moment, d'éviter la saturation des réseaux et d'anticiper en cas 
de fortes pluies évitant ainsi les débordements fortement polluants. Il offre une vision 
globale du fonctionnement de l’ensemble du système d'assainissement à partir de 
données métrologiques et météorologiques. Cette solution est capable de réaliser un 
scénario tendanciel d'évolution hydraulique du système d'assainissement, dans les 


24h à venir par temps sec, ou dans les 2 à 6h à venir par temps de pluie. 


Bref, AQUADVANCED® met en œuvre une modélisation hydrologique et 
hydraulique des ouvrages, alimentée en temps réel par des données de temps sec 
et des données de pluies prévues et mesurées par radar météorologique. La 
simulation hydraulique est couplée à un modèle de dispersion de la pollution 


rejetée en mer par les déversoirs d'orage. Ce couplage permet de proposer les 
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stratégies de gestion et les scénarios de pilotage les mieux adaptés au contexte 
grâce à une vision prédictive de la pluie et de la sensibilité des milieux 
récepteurs. AQUADVANCEDS fournit ainsi une aide à la protection du milieu 
récepteur, à la gestion active des plages, mais aussi à la gestion du risque 


d'inondation. C'est une technologie innovante basée sur le contrôle optimal . 





IX. Conclusion 


La « gestion de référence » d'un réseau d'assainissement permet de renseigner 
l'exploitant sur l'opportunité de gérer un réseau d'assainissement en temps réel 
dans l'objectif, par exemple, de minimiser les volumes déversés. Parmi les 
techniques de recherche opérationnelle, la programmation linéaire par les 
graphes et la programmation linéaire mixte permettent de fournir les stratégies 
de commande des organes de contrôle (pompes, vannes motorisées). Nous 
avons montré qu'il est possible de modéliser de façon simplifiée par les deux 
techniques précédentes les ouvrages disponibles en assainissement tels que les 
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collecteurs, bassins et déversoirs d'orage. Toutefois, il est nécessaire 
d'introduire des variables binaires dans le cas de la programmation linéaire 
mixte et un ordre de saturation des arcs par l'approche par les graphes. 
L'algorithme utilisé par la programmation linéaire mixte « Branch and Bound » 
peut fournir une solution minimale globale moyennant un temps de calcul qui 
peut être très important. Par contre, bien que pour la programmation linéaire 
par les graphes le temps de convergence soit très rapide, on ne garantit plus un 
minimum global. En effet, l'algorithme utilisé impose une indépendance sur le 
choix des arcs à saturer, ce qui n'est pas vérifié. Compte tenu des avantages et 
inconvénients de chaque approche, nous avons voulu quantifier sur l'exemple 
du réseau d'assainissement de Saverne le temps de calcul et les différences en 
terme de volume déversé. Les résultats ont montré que la programmation 
linéaire mixte fournit des temps de calcul qui peuvent atteindre plus de 24 
heures. Par contre, l'approche par les graphes permet un temps de calcul de 
l'ordre de 5 minutes en moyenne avec un minimum global en terme de volume 
déversé atteint qui n'excède pas 5 % par rapport à la solution fournie par la 
programmation linéaire mixte. On vient ainsi de montrer qu'on peut choisir un 
compromis entre la détermination exacte d'une optimisation et le temps de 
calcul qui lui est associé avec une solution approchée acceptable présentant un 


temps de calcul beaucoup plus faible. 
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